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ABSTRACT 

A  simple  analytic  model    is   formulated  which   should   prove  useful    for 
evaluating   a  satellite  surveillance  system.      This  model   describes   the 
fraction  of  orbits   in  which  a  long   lived   satellite  will    "see"  a   selected 
position  on  the  earth's   surface.      The   fraction  of  orbits,   denoted   by  P, 
is   a  function  of  the   latitude  of  the   ground   position   in  question,    the 
orbital    inclination  of  the   satellite,    and  the   sweep  width  of  the  satel- 
lite's  sensor. 

P  is  used  to  find  estimates  for  the  mean  time  to  next  sighting  of  a 
selected  ground  position  given  that  position  has  suddenly  become  of  in- 
terest, to  estimate  the  number  and  locations  of  ground  readout  stations, 
and  to  estimate  the  number  of  times  a  satellite  would  sight  a  transitting 
object  on  the  earth's  surface.  In  addition,  the  application  of  P  to  the 
problems  of  optimum  orbit,  evasive  counter-measures ,  and  satellite  weapon 
delivery   is   discussed. 
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1.   Introduction. 

Rapid  developments  in  space  technology  have  resulted  in. an  increased 
use  of  unmanned  earth-orbitting  satellites  for  both  civilian  and  military 
purposes.   These  satellites  provide  information  on  outer  space,  the  earth's 
atmosphere,  and  the  earth's  surface.   Satellites  for  earth  surveillence 
offer  a  distinct  advantage  over  earth-bound  systems  because  of  their  abil- 
ity to  scan  large  segments  of  the  earth  in  a  relatively  short  period  of 
time,  and  because  they  are  unencumbered  by  national  air-space  restrictions. 

The  Tiros  and  Nimbus  series  weather  satellites  represent  well  known 
examples  of  the  use  of  satellites  by  the  United  States  for  earth  surveil- 
lence.  Tiros  III,  launched  in  1961  was  equipped  with  photographic  sensors 
and  provided  information  on  storms  and  storm  paths.   Nimbus  I,  launched  in 
1965  served  a  similar  purpose  but  made  use  of  television  cameras  as  sen- 
sors. 

Although  weather  received  the  first  attention,  it  is  not  the  only 
phenomenon  worthy  of  satellite  surveillence.   Nuclear  test  explosions  are 
of  sufficient  interest  that  two  Vela  Hotel  satellites  were  launched  by 
the  United  States  in  1965.   Additional  satellites  apparently  have  been 
launched  for  certain  classified  military  purposes.   The  Samos  series 
satellites  evidently  have,  for  the  last  five  years,  been  photographing 
military  installations  in  Russia  and  Red  China.  [6j  Seven  military  recon- 
naissance satellites  (possibly  of  the  Samos  series)  were  reported  to  have 
been  launched  during  1965.  [V|   These  satellites  had  polar  orbits  and 
mission  lives  of  about  five  days.   The  observations  made  by  these  satel- 
lites are  unknown,  but  could  include  ship  and  troop  movements  or  test 
firings  of  ICBM's,  as  well  as  military  installations. 

Given  a  surveillence  task,  alternative  satellite  systems  might  be 
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designed  to  perform  this  task.  We  may  want  to  decide  which  system  or  com- 
bination of  systems  is  most  cost-effective.   Or  we  may  want  to  compare 
the  alternatives  with  more  conventional  systems  such  as  U2  flights.   To 
make  such  comparisons  we  must  be  able  to  evaluate  the  effectiveness  of 
each  type  of  surveillence  system.   This  paper  will  be  concerned  with  eval- 
uating satellite  surveillence  systems. 

What  elements  are  important  in  the  evaluation  of  a  satellite  surveil- 
lence system?   This  is  a  difficult  question  to  answer.   The  appropriate 
elements  must  include  both  the  nature  of  the  phenomena  to  be  observed,  and 
the  available  hardware,  including  the  satellite,  its  components,  and  ground 
support  equipment.   The  nature  of  the  phenomena  to  be  observed  may  be  such 
that  one  satellite  of  a  particular  type  is  sufficient.   Or,  it  may  be 
found  that  multiple  satellites  of  one  or  more  types  are  required.   The 
appropriate  number  will  depend  upon  the  effectiveness  of  each  individual 
satellite.   Thus,  an  initial  appraisal  of  a  system  should  be  based  on  an 
individual  satellite's  performance. 

What  features  of  an  individual  satellite's  performance  should  be  ex- 
amined?  Two  features  appear  to  be  dominant.   For  those  phenomena  which 
occur  randomly  in  time  and  location  on  the  surface  of  the  earth,  it  seems 
desirable  to  make  some  statement  about  the  expected  proportion  of  the  total 
phenomena  of  interest  the  satellite  will  observe.   For  other  phenomena  we 
are  more  interested  in  seeing  the  result  of  an  event  rather  than  the  event 
itself.   Thus,  we  would  like  to  know  the  average  elapsed  time  before  a 
satellite  observes  a  location  which  has  suddenly  become  of  interest.   (For 
example,  if  on  the  average  100  natural  disasters  occur  in  a  year  in  the 
area  covered  by  the  satellite,  then  we  may  want  to  know  the  average  time 
lapse  between  occurrence  of  a  disaster  and  observance  of  the  place  of 


disaster.) 

The  present  procedure  for  satellite  system  evaluation  is  to  construct 
an  elaborate  model  describing  the  location  of  the  satellite  at  all  times 
for  a  given  orbital  inclination,  altitude,  etc.   This  approach  requires 
that  every  combination  of  orbital  inclinations,  sweep  width,  altitude, 
etc.,  must  be  investigated  to  have  a  complete  data  bank  for  estimating 
the  frequency  of  observation  or  the  time  before  an  event  is  observed. 
This,  by  nature  of  the  model's  complexity,  requires  the  use  of  a  high 
speed  computer  and  much  computation  time.   For  example,  if  each  of  three 
parameters  is  to  be  investigated  for  10  different  values,  then  1000  separ- 
ate replications  would  be  required. 

The  author  believes  that  an  alternative  approach,  by  using  a  simpli- 
fied analytical  model,  may  prove  to  be  just  as  useful  for  most  appraisals 
of  satellite  surveillence  systems.   Appropriate  analytical  expressions, 
even  if  approximations  to  the  actual  situation,  usually  help  to  point  out 
the  important  parameters,  variables,  etc.  As  a  consequence,  a  simplified 
model  like  that  discussed  in  this  paper  may  point  out  some  possible  con- 
clusions not  immediately  apparent  in  a  computerized  model.   For  example, 
it  was  mentioned  that  the  United  States  launched  seven  surveillence  satel- 
lites during  1965  which  were  characterized  by  polar  orbits.   It  will  be 
obvious  from  the  model  to  be  presented  in  this  paper  that  for  reconnaissance 
purposes  a  polar  orbit  is  extremely  inefficient  unless  the  polar  regions 
themselves  are  of  particular  interest. 

Basic  to  the  approach  is  the  development  of  a  fairly  simple  analy- 
tical expression  for  the  likelihood  of  seeing  any  event  of  interest.   It 
would  be  convenient  if  we  could  find  the  probability  that  the  satellite 
passes  over  a  chosen  ground  position  during  any  randomly  selected  orbit— 


but   since  the  satellite's   position  is  completely  deterministic  we   can  not 
speak  of  the   "probability"  that  the   satellite  passes  over  such  a  position. 
(For  a  selected  orbit  either   the  satellite  passes  over  the  position  or   it 
doesn't.)     However   it  does  make  sense  to   speak  of  the   "fraction     of  orbits" 
in  which  the   satellite  passes   over  a   selected  ground   position.      This  dis- 
tinction between  "fraction"  and  "probability"   is   a   fine  one,    but   an   impor- 
tant one.      For  example,    if   a   satellite  sees  London   every  seventh  orbit   then 
it   sees   London  on  one-seventh  of   the  orbits,    but  the  probability   it   sees 
London  on  any  randomly  selected  orbit    is   not  one-seventh. 

The   likelihood  function  to  be  derived  will,    then,   describe  the  frac- 
tion of  orbits   that  the  satellite  spends    in  "positive"  orbits,    where  a 
positive  orbit   is  defined  to   be  one  such  that  the  satellite  passes  over 
the  ground  position  in  question  on  that  orbit.      If,    for  example,    the  ground 
position   is  Beirut,   then  the  fraction  of   positive  orbits  would  be  the  ratio 
of  the  number  of  orbits    in  which  Beirut    is  within  the   sensors   sweep  width 
to  the  total   number  of  orbits    (per  some  time  interval,    e.g.    the   life  of 
the   satellite).      If  P  represents  this   fraction,   then,    using  reasonable 
assumptions,    P  can  be  shown  to   be  a   function  of  the  satellite's  orbital 
inclination,   the  sweep  width  of  the  sensor,    and  the  latitude  of   the  point 
of   interest.      P  is   also  a  function  of  the  satellite's   altitude,    in  that 
the  sweep  width   is   usually  a  function  of  the   satellite  altitude. 


2.  Formulation  of  the  Likelihood  Model 

The  formulation  of  the  likelihood  model  is  facilitated  by  the 
following  sketch. 


^  Satellite 


Satellite  motion  in  general  is  confined  to  being  great  circle  mo- 
tion.  That  is,  the  ground  track  of  the  satellite  describes  a  great 
circle  on  a  stationary  globe.   Controllable  parameters  are  altitude, 
orbital  inclination,  and  sweep  width. 

For  surveillence  purposes,  upper  and  lower  bounds  exist  on  satel- 
lite altitude.   The  lower  bound  is  caused  by  the  earth's  atmosphere.   A 


satellite  must  be  outside  the  atmosphere  to  remain  in  orbit.   The  upper 
limit  is  dependent  on  the  resolution  capabilities  of  the  satellite  sen- 
sors.  The  lower  and  upper  bounds  are  presently  about  100  and  300  miles 
respectively.   The  300  mile  maximum  altitude  means  that  the  time  for  one 
orbit  will  be  less  than  1.4  hours.   That  is,  while  the  earth  rotates  one 
revolution  the  satellite  has  orbitted  about  17  times.   This  formulation 
will  assume  that  altitude  and  sweep  width  are  directly  related  and  that 
knowing  sweep  width  is  equivalent  to  knowing  satellite  altitude.   There- 
fore no  further  mention  of  altitude  need  be  made. 

The  orbital  inclination  angle  -Q-   is  equal  to  the  extreme  latitude 
attained  by  the  satellite.   Once  launched,  we  assume  *0"  remains  unchanged 
for  the  life  of  the  satellite. 

The  sweep  width  2r  is  the  width  of  the  path  on  the  earth  which  the 
satellite  sensor  sweeps  out  as  it  proceeds  on  its  orbit.   For  photographic 
sensors  2r  would  be  the  width  of  the  area  photographed  on  one  exposure, 
for  radar  it  would  be  the  length  of  one  scan,  scans  being  perpendicular 
to  the  direction  of  the  satellite  motion. 

We  let  o<  represent  the  latitude  at  which  an  event  occurs,  and  since 
we  have  no  control  over  the  location  (or  time)  of  events  of  interest,  <* 
is  an  uncontrolled  parameter.   The  ground  position  at  latitude  of.    which 
is  of  interest  will  be  designated  as  A,  where  A  is  either  selected  ran- 
domly or  non-randomly,  depending  on  the  context. 

General  Assumptions  of  the  Model 

ITie  following  general  assumptions  are  basic  to  the  likelihood  model: 

(i)    The  satellite  under  consideration  remains  in  orbit  for  long 
periods  of  time,  upwards  of  1000  orbits.   The  life  of  a  sat- 
ellite will  be  denoted  by  L. 
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(ii)   An  event  of  interest  is  equally  likely  to  occur  at  A  at  any 
time.   Given  an  event  occurs  during  the  time  L,  then  we 
assume  the  random  variable  denoting  the  time  of  occurrence 
is  uniform  on  (0,L). 

(iii)   0  cannot  be  varied  after  launch. 

(iv)   The  orbit  is  circular. 

(v)  Every  point  on  earth  which  lies  within  the   sweep  width  of  the 

satellite's   sensor  as  the  satellite  passes  overhead   is  seen 
with  probability   1. 

(vi)        The  earth  is  a  perfect   sphere. 

(vii)      The  sensor  points  normal   to  the  earth's  surface  at  all   times, 
and  can  see  a  distance  r  to  each  side  of  the  satellite's 
ground  track.      This  varies  with  the  actual   sensor  aspect   in 
designed  systems    (in  general   the  sensors  are  aimed  to  the 
side)    but   simplifies  the  trigonometry,   without  changing  any 
basic  concepts  of  the  system. 

(viii)    Sensors  are  either  photographic  or  radar.      In  either  case, 
the  ensuing  model   assumes   perfect  visibility.      (Brief  men- 
tion is  made  in  Section  4  of  how  one  might   incorporate  the 
weather  problem,    and  consequently  the  problem  of  general 
atmospheric   attenuation,    into  the  modeL) 
The  Question  of  Accuracy 

The  formulation  to  be  presented  does  not  yield  extreme  accuracy  of 
results.      Extreme  accuracy  seems  unnecessary  in  the  problem  of   initial 
system  evaluation.      For   example,    if  we  determine,   through  the  model,   that 
the  mean  time  to   next  sighting  of  a  ground  position  A  is   2  hours  when  in 
actual   fact   it   is  2.078  hours,   then  this   probably  would  be  of  small  con- 
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sequence  in  any  decision  made  regarding  the  system.   Certainly  it  would 
be  misleading  to  present  results  implying  accurately  to  six  or  seven  fig- 
ures when  basic  assumptions  preclude  accuracy  to  more  than  two  or  three 
figures.   Computer  formulations  specifically  designed  for  the  surveillence 
evaluation  problem  may  claim  accuracy  to  six  figures.  [33  These  figures 
may  or  may  not  be  that  accurate  when  compared  to  an  actual  satellite  in 
orbit,  but  any  initial  evaluation  of  a  satellite  system  would  probably  have 
remained  the  same  had  the  figures  been  rounded  to  two  significant  figures, 
which  is  approximately  the  accuracy  inherent  in  the  proposed  model. 

The  Likelihood  Function 

We  will  denote  the  "fraction  of  positive  orbits"  by  P(#  ,  <*  ,  r) . 
However  V(&  ,  o(  ,  r)  will  often  be  referred  to  as  the  likelihood  function 
P.   As  the  notation  implies,  P  depends  on  three  parameters:  orbital  in- 
clination 9,  latitude  o<  of  the  ground  position  A  to  be  observed,  and  the 
sweep  width  2r. 

The  analytical  expressions  for  P(£  ,  o<  ,  r)  are  derived  in  Appendix 

I  and  are  given  below.   The  derivation  is  based  on  the  assumption  that  P 

is  equivalent  to  the  ratio  of  the  sweep  width  measured  along  a  latitude 

to  one-half  the  circumference  of  the  latitude.   (i.e.  TS   in  the  previous 

LL* 
sketch) 


P(*.  «;r)  =  \  /-  Jr  {o>*-'[~( s*  i 


+■  s/no(  cos 


*s 
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#-g  £  «  <  e+S 
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where  S  =  r  and  R  is  the  radius  of  the  earth. 
R 

Figure  1  illustrates  the  behavior  of  P  as  a  function  of  o<  for 
several  values  of  6  when  r  =  150  nautical  miles.   The  figure  shows 
that  P  attains  its  maximum  when  o(  =  &- S  .     This  fact  can  be  con- 
firmed by  an  investigation  of  (1). 
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3.      Applications. 

Mean  of  Time  Until  Next   Sighting 
One  of  the  present  parameters  used  for  appraising  a  satellite  sur- 
veillence  system  is  the  mean  time  between  sightings.1     From  our  model  the 


TV 


mean  time  between  sightings  is  simply  _-£-  -   However  if  location  A  has 
suddenly  become  of  interest  we  probably  would  not  care  about  the  time 
lapse  between  the  last  sighting  of  A  and  the  next  one,  but  rather  on  the 
time  lapse  from  the  present  until  the  next  sighting.   These  two  parameters 
are  by  no  means  equivalent. 

In  this  section  we  will  develope  several  estimates  for  the  mean  time 
from  a  given  instant  until  the  next  sighting  of  a  specified  point  A.   In 
these  developments,  we  will  assume  that  phenomena  of  interest  are  equally 
likely  to  occur  at  A  at  any  time  during  the  life  of  the  satellite. 

If  T  is  the  time  lag  between  occurrence  of  an  event  at  location  A 

Li 

and  the  first  sighting  of  A  after  the  occurrence,  then  we  are  looking  for 
the  mean  value  of  TL,  denoted  by  M(TL) .  An  exact  expression  for  M(TL>  is 
difficult  to  obtain,   but  various  estimates   can  be  made. 

To  see  why  the  exact  expression  is  difficult  to  obtain  analytically 
consider  the  following  sketch. 


'■For  an  example,    see   03- 
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If  we  consider  an  upward  (downward)  pass  through  the  latitude  o(  to  be  one 

such  that  the  satellite  is  moving  from  o<  to  a  more  northerly  (southerly) 

latitude,  then  the  satellite  can  see  point  A  on  any  given  orbit  in  one  of 

two  ways.  Either  it  sees  A  on  an  upward  pass  or  on  a  downward  pass.   The 

time  between  an  upward  and  downward  pass  depends  on  the  latitude  o<  .   If 

T*  deontes  this  time  and  T  is  the  time  for  a  satellite  to  make  one  com- 

o 

plete  orbit,  then  the  discrete  points  in  time  at  which  sightings  can  pos- 
sibly be  made  occur  in  a  sequence  T,  ,  T2, ,  T., ,  where  the  times 

between  sightings  are  t' ,  TQ  -  T* ,  t' ,  TQ  -  T* 

In  reality,  sightings  do  not  occur  in  a  neat  sequence  T, ,  T~,...,  T^, 
...  as  described  above.   That  sequence  merely  tells  when  sightings  are  pos- 
sible.  The  actual  sightings  may  occur  in  sequences  which  are  exceedingly 
difficult  to  describe  analytically.   For  example,  as  indicated  by  the 
previous  sketch,  a  sighting  occurs  when  on  an  upward  pass  the  satellite 
passes  over  point  P  while  A  lies  between  Q  and  Q  .   Since  the  distance 
QQ'  is  on  the  order  of  several  hundred  miles,  the  speed  of  the  earth's 
rotation  makes  it  possible  for  A  to  transit  QQ*    in  much  less  time  than  it 
takes  the  satellite  to  make  a  complete  orbit.   Thus  the  satellite  may  have 
just  crossed  o<  while  A  was  just  west  of  Q,  and  A  will  pass  Q*  before  the 
satellite  makes  another  upward  pass.   A  "sighting  possibility"  would  thus 
be  missed. 

Rather  than  attempting  the  difficult  task  of  predicting  the  actual 
sequence  of  positive  orbits,  we  will  develope  estimates  for  M(T-r)  .   If 
tests  show  these  estimates  to  be  reasonably  accurate,  then  they  should  be 
almost  as  useful  to  the  decision  maker  as  the  true  value.   Let  us  consider 
three  possible  estimates. 
Estimate  I:   If  we  assume  actual  sightings  occur  at  equally  spaced  inter- 
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vals  of  time,  then  the  spacing  can  be  obtained  from  the  likelihood  function 
P.   The  assumption  that  events  of  interest  occur  with  equal  likelihood  at 
any  time  during  the  life  of  the  satellite  implies  that  if  an  event  occurs 
during  the  time  (T^,  Ti+1)  then  it  is  equally  likely  to  occur  at  any  in- 


stant in  that  interval;  that  is,  the  conditional  time  of  occurrence  T  is 


T 
uniform  on  (T^  Ti+1>.   Since  the  length  of  this  interval  is  — 2—  ,  M(T  ) , 


an  estimate  for  M(T-r )  ,  is 

T 
M(T,  )  =  — 2-  . 
L     2P 

Intuition  suggests  that  this  may  be  a  good  estimate  near  the  equator 
where  the  time  between  an  upward  and  downward  pass,  and  therefore  the  time 
between  possible  sightings,  is  the  constant  %T  .   Since  the  times  of  pos- 
sible sightings  are  equally  spaced,  the  time  between  actual  sightings  may 
tend  to  be  equally  spaced  also.   But  it  seems  equally  intiutive  that  as  A 
moves  away  from  the  equator  the  estimate  will  become  less  reasonable. 
Estimate  2;   Arbitrarily  define  a  probability  model  for  the  orbiting  sat- 
ellite, wherein  P,  instead  of  being  the  fraction  of  positive  orbits,  is 
defined  to  be  the  probability  of  a  positive  orbit.  Further  assume  that 
all  orbits  are  independent  of  each  other.   This  model  is  then  equivalent 
to  one  with  Bernoulli  trials  with  probability  of  success  equal  to  P.  Fur- 
thermore the  number  of  trials  to  first  success  for  such  Bernoulli  trials 
is  a  geometric  random  variable  with  mean  ——,  where  P  is  the  probability  of 
success  on  each  trial.   Equivalently,  in  terms  of  our  model,  the  number  of 

orbits  to  the  first  success  (positive  orbit)  is  — L_  and  therefore  the 

T 
estimate  for  the  time  to  first  positive  orbit  would  be      . 

Intuition  suggests  that  this  estimate  should  be  good  when  o(  ap- 
proaches 9,  but  not  so  good  near  the  equator.  The  time  between  an  upward 
and  downward  pass  approaches  zero  as  oC-^-0*,  and  between  downward  and  up- 
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ward  approaches  T  .   Thus  the  time  till  possible  sighting  will  be  twice  as 
long  as  that  at  the  equator. 

Estimate  3:   Since  neither  of  the  above  estimates  seems  reasonable  for  all 
values  of  o(  >  we  return  to  the  physics  involved  to  see  if  an  assumption 
can  be  made  which  provides  abetter  estimate. 

If  a  ground  position  A  is  at  latitude  &    we  have  said  it  can  be  seen 
during  that  orbit  at  one  of  two  times,  it  can  be  seen  during  the  upward 
pass  or  the  downward  pass.   It  is  relatively  easy  to  obtain  an  exact  ex- 
pression for  T  ,  the  time  between  the  upward  and  downward  passes.   Knowing 
T  and  T  we  can  describe  the  sequence  of  times  at  which  sightings  are 
possible.   We  also  know  that  the  fraction  of  positive  orbits  is  P. 

Using  these  two  pieces  of  information  we  assume  that  actual  sightings 
occur  in  a  sequence  such  that  the  time  between  sightings  is  proportional 
to  the  time  between  possible  sightings,  where  the  constant  of  proportion- 
ality is  such  that  the  fraction  of  positive  orbits  is  P.   We  get  (see 
Appendix  II) 


Ami*  £p  [[cos-'fco/  &  tan  «)f  + 


[rr  -  cos'Yco/-  &  fan  «$f] %        (  3 ) 


Notice  that  for  o(  =  0,M(TT)  becomes  — _  which  is  identical  to 

L  2P 

A  T 

estimate   1,    and   for    <X    =   9 ,   M(TT  )   becomes        °       which    is    identical   to 

l  p 

estimate  2.  Thus  estimate  3  varies  with  o<  and  also  fits  what  intuition 
suggests  should  be  the  case  at  the  extremes  of  the  range  of  o(  .  Further 
evidence  of  this  estimate's   reasonability   is  given  by  Appendix  III  where 

A 

M(Tt)  is  compared  with  results  obtained  from  a  computerized  simulation 


developed  by  the  Naval  Missile  Center.  £3] 
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Location  and  Number  of  Readout 
Stations  Required  for  a  Specified  Level  of  Contact 

Any  study  of  the  effectiveness  of  a  satellite  surveillence  system 
must  include  an  appraisal  of  problems  of  transferring  'the  information 
from  satellite  to  evaluation  centers  located  on  earth.   To  consider  a  sat- 
ellite which  orbits  for  long  periods  of  time  we  should  discuss  readout 
stations.   Readout  stations  are  radio  receiving  stations  capable  of  re- 
ceiving signals  from  a  satellite  whenever  the  satellite  is  within  a  cer- 
tain distance  of  the  readout  station.  We  will  refer  to  the  circular  area 
wherein  such  reception  is  possible  as  the  transmission  area.  As  the  sat- 
ellite passes  through  a  transmissidn  area  it  transmits  information  accumu- 
lated since  the  last  pass  through  a  transmission  area.   The  author  knows 
of  no  present  study  which  analyzes  the  optimum  location  of  such  stations 
or  the  number  needed  to  ensure  adequate  transmission  time,  therefore  it 
seems  worthwhile  to  consider  this  question. 

P(£,  o(  ,   r)  can  be  used  to  determine  the  percentage  of  time  a  satel- 
lite will  be  in  contact  with  a  particular  readout  station.   The  procedure 
we  will  use  is : 

1.  Find  the  diameter  of  the  transmission  circle  within  which 
radio  signals  from  the  satellite  are  receivable. 

2.  Estimate  M(t) ,  the  mean  time  a  satellite  takes  to  cross  the 
transmission  circle  given  it  intersects  the  circle. 

3.  Determine  the  value  P(£-,  crt  ,  r)  ,  where  o(  is  the  latitude 
of  the  transmission  circle  and  r  is  the  radius  of  the  transmission  circle. 

4.  Estimate  the  mean  time  per  orbit  available  for  transmitting, 
by  P  %  M(t)  . 

Knowing  the  rate  of  transmission  of  information,  a  reasonable  esti- 
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mate   for  the  number  of   readout   stations   needed  to  maintain  a   specified 
level   of   information  transmission  can  be  obtained. 

To   find  the  diameter  of   the  transmission  circle  consider  the   fol- 
lowing sketch. 


In  this   sketch  h   is  the  height  of  the   satellite  above  the  earth's   sur- 
face,   and   a   is  the  minimum  reception  angle  above  the  horizon. 
From  the   law  of   sines, 

/?  =  R+H  R+H 

-*  //?  cos  a  ) 


There&re^       AnjkABC=      &$'[[££ 


R>+  H 


an 


</ 


t 


cos 


-l  /  R  cos  a 


-  a  . 


Thus,  the  diameter  of  the  transmission  circle,  d,  is 


c/=  Z(#+V/*    . 


w 
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Thus  we  see  that  a  readout  station  can  receive  signals  from  the  satel- 
lite if  the  satellite  is  within  a  circle  with  diameter  2(R  +  h)0  with 
center  at  the  readout  station.   If  we  nov  let  r  in  the  expression  for 
P(-&,  <*  ,  r)  equal  1(R  +  h)0  ,  then  we  will  have  the  expression  for  the 
fraction  of  orbits  wherein  the  satellite  passes  within  the  transmission 
circle. 

To  find  the  mean  time  available  for  transmitting  given  the  satellite 
passes  through  a  transmission  area,  we  will  assume  that  the  distribution 
of  paths  which  the  satellite  makes  as  it  crosses  the  transmission  area  is 
uniform;  that  is,  the  satellite  is  equally  likely  to  cross  the  circle  on 
any  chord.   Finding  the  mean  length  of  path  based  on  this  assumption  is  not 
difficult.   If  the  radius  of  the  transmission  circle  is  r  then  the  mean 
distance  travelled  while  in  the  circle  is  approximately  ^r   •  (See  Appen- 
dix IV). 

From  the  above  results,  we  can  obtain  an  estimate  for  the  mean  time 
per  orbit  spent  in  the  transmission  area  of  a  single  readout  station.  It 
is 

M(t)  =  -^-  P(e,  <*  ,  r);    r  =  (R  +  H) 

V  =  ground  velocity  of  satellite. 

Because  P  is  maximized  when  o(  -   &  -  £   and  M(t)  is  linear  in  P,  we 
can  maximize  the  mean  transmission  time  per  orbit  by  placing  the  readout 
station  at  latitude  -0-  -  &  . 

We  next  consider  the  question  of  how  many  readout  stations  are  re- 
quired to  attain  at  least  some  specified  mean  transmission  time  per  orbit. 
We  can  reasonably  assume  that  all  stations  would  be  placed  close  to  lat- 
itude &  -  h    •   We  also  assume  the  stations  are  spread  out  on  that  lati- 
tude to  avoid  a  satellite  passing  through  two  stations  in  rapid  succession 
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and  then  not  passing  through  another  for  many  orbits.   Under  these  assump- 
tions a  reasonable  estimate  for  the  number  of  readout  stations  N  required 

to  average  t'  minutes  of  transmission  per  orbit  is 

i 

N  =  j.\    g   • 
B(t) 

If   a  value  of   t      is   designated   as   the  minimum  mean  time  per  orbit   needed 
in  order  to   receive  some   specified   fraction  of   the  total    information  ac- 
quired  by  the   satellite,    then  N  tells   us   how  many   readout    stations   are 
required.      Such   information  would  be   needed   to   estimate  the  cost  of   a   pro- 
posed  satellite  system. 

Estimating  the  Number  of   Intersections  of   a   Trans itting  Ship 

by  a  Satellite 

To   find   an  estimate   for   the   number  of    intersections  of   a  moving  point 
on  earth,    consider  a  ship  going  from  point   a  to   destination  b,    where  a 
and  b  are  at    latitudes    A    and  o    respectively.      Let   us   assume: 

1.  The  ship  proceeds   at   constant   speed   s. 

2.  The  ship  moves  on  great   circle  path. 

3.  The  distance  d   from  a  to  b  on  this   path   is   known. 

4.  The  track  of  the  ship  does   not   leave  the   area  of   satellite 
coverage. 

First,   to   avoid   complicated  expressions  derived   from  differential   equa- 
tions  giving  the  instantaneous    latitude  of  the  ship  as    it   proceeds,   we 
elect  to   use   a  constant,   }f   ,    to   represent  the  time   averaged  mean   latitude 
of  the  ship;    that   is ,   2T     is   an  estimate  for  the   latitude  which,    if  used  to 
replace   o<    in  P(&,o(    ,    r)  ,   would  give  the   same  result   as   if   we  continu- 
ously evaluated  P  at   each   instant  of   time  during   the  ship's   transit,     tf    is 
defined  as 
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Since  the   ship's  track   is  a  great  circle  track  and  the  ship's  speed   is 
constant,    this   is  a  lower  bound  for  the  ship's   time  averaged  mean  latitude. 
A  conservative  estimate  for  the  number  of   intersections   will   consequently 
be  obtained. 

The  derivation  of  the  estimate  of  the  total   number  of  sightings   for 
a  ship  trans itting  from  a  to  b  under  the  above  assumptions   is  the  follow- 
ing: 

Let     h  =  satellite  altitude 

t   =   ship's  transit  time  from  a  to  b 

v  =  satellite's   speed 

g   =  gravitational   constant. 


Then 


Now 


*=/ 


£ ,   and  t   =  -1 


(R  +  h) 


T     -     2TT(R  +  h) 
o   "  v 

and  thus  N,   the  number  of  orbits  during  the  transit   is 
t 


N   = 


To 


If   I   is  the  number  of    intersections  during  the  ship's   transit   then  an 
estimate  for   I   is 

I  =  NPC-0-,  X  ,   r). 
Example:      If  the  ship  goes   1000  miles   from  10°   to   20°   at   12   knots,   the 
average  number  of  sightings  by  a   200  mile  high  satellite  with  orbit 
6  =   50°   and   r  =   150  will   be 

I    =   P    •    N  =  NP(£,ff    ,    r)    =    (63.2)(.035)    =   2.2. 
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4.   Areas  for  Further  Study. 

Having  derived  an  expression  for  the  fraction  of  positive  orbits  one 
can  investigate  various  areas  other  than  those  already  mentioned.   The 
following  are  areas  which  seem  to  be  of  obvious  interest  in  the  evaluation 
task. 

1.   Finding  the  optimum  orbit  based  on  some  measure  of  effectiveness 
is  a  problem  of  considerable  interest.   Figure  2  shows  graphs  indicating 
the  distribution  of  major  nautical  activities  as  a  function  of  latitude. 
Further  work  could  be  accomplished  correlating  P  and  these  distributions 
to  optimize  the  orbit.   For  example,  if  the  measure  of  effectiveness  is 
"number  of  contacts  per  unit  time"  then  we  could  proceed  as  follows.   Des- 
cribe a  graph  of  figure  5  by  a  function  f(o<),  where  f  is  either  a  discrete 
or  continuous  function  approximating  the  graph.   If  we  consider  southern 
latitudes  to  be  negative,  then  we  can  let 


F(«)  =  /  f(«')d*'        or       FM  =H  *(«')  , 


depending  on  which  is  applicable.   Then,  let 


G(&L>r)   =j'pf*>«)r)c/f(«).  (io) 


If  G   is  the  maximum  of  G  with  respect   to   0^,    then  that  9^  which  yields 
G   is   the  optimal  orbital    inclination  for  our  chosen  measure  of  effective- 
ness. 

The  distribution  of  naval   activities   with  longitude  may  be   combined 
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Figure  2.  Distribution  of  Various  Naval  and  Maritime  Activities 
with  Respect  to  Latitude.  [2J 
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with  the  distribution  for  latitude.   The  optimum  orbit  for  observing  speci- 
fied regions  may  then  be  quite  different  than  if  latitude  distribution 
only  are  considered. 

2.  In  studying  the  effectiveness  of  a  satellite  as  a  weapon  delivery 
vehicle,  we  need  only  replace  the  half  sweep  width  r  by  the  lethal  range 
r'  of  the  weapon  carried  by  the  satellite.   Then,  if  the  ground  position 
of  the  satellite  is  within  the  lethal  range  r'  of  the  target,  the  weapon 

A 

could   be  dropped.      P  and  M(T-, )    would  have  the   same   form  as   above,    but 
M(TL)    would  be  the  mean  time  until    a  weapon  could  be  dropped  on  a   target 
after  the  decision  to   destroy  the  target   had   been  made.      One  might   cor- 
relate this  with  the   possibility  of    inflight   satellite  course  changes. 

3.  Investigating  counter-measure  techniques  may  prove  of    interest.      In 
this   respect   figure   1   gives   some    indication  of   P  and   its    sensitivity  to 
latitude   in  the  extreme   regions   near  9.      Knowing  the   inclination  of   an 
enemy's  orbiting  satellite  and  correlating  that    information  with  the   impli- 
cation of   figure  1  would   indicate  ship  routing  procedures   we  could   use  to 
avoid  detection.      We  obviously  would  try  to   route  our  ships   to   stay  clear 
of  the   latitudes   near     ■&  -  &     ,    where   P  is   largest. 

4.  When  data   is   available  on  global   cloud   cover  distributions,    P  can 
be  modified   by  the  probability  of   cloud  cover   and/or  darkness  obscuring 
the   sensor's  vision   (as    in  the  case  of  photographic  sensors).      As   a  con- 
sequence,   P  would  become  closer  to   being   a  probability  than   a  fraction. 
As   a  result  of  cloud  cover  being   a  function  of   longitude  as   well   as   lat- 
itude,   P  would  also   have  to  be  dependent   on   longitude.      The  problem  of 
predicting  cloud  cover    is  obviously  complex  and   any  model   developed   to 
describe  the  distribution  of  cloud  cover  would  probably  have  to   be  a  com- 
puter simulation. 
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APPENDIX   I 
DERIVATION  OF  THE   FRACTION  OF  POSITIVE  ORBITS 

The   function  ?(&  ,  <  ,    r)    is  defined  to   be   the  fraction  of 
positive  orbits  to  total  orbits  given  a  satellite   in  orbital    inclina- 
tion 9,    a  sensor  sweep  width  2r,    and   a  ground   position  at    latitude  c<  . 
The  basic  rationale   for  finding  P   is   to   consider  the  circle  on  the 
sphere  made  by  latitude  <tf  ,    and  to   find   the  ratio  of  the  path  width  of 
the  satellite  measured  along   latitude   u    to  one-half   the  total   circle. 
Then  if  the  ground  position  A   is   on   latitude  c*  ,    the  fraction  of  times 
A  is   seen  will   be  proportional   to   this   ratio. 

To  obtain  an  analytic   expression  for   P  based  on  this   rationale 
consider  the  following   sketch,    in  which  capital    letters   denote  points 
on  the  sphere,    lower  case   letters  denote  angles  of   spherical    triangles, 
and  Greek   letters  denote  sides   of  spherical   triangles,    a  side  being  the 
usual   radian  measure  of  the  central   angle  subtended  by  that   side. 

N 
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If  we  could  determine  the  length  of  BB',   then  divide  BB'    by  FG,    (one 
half  the  circumference  of  the  circle  made  by  latitude  t*  )   we  would  have  P. 
However   it    is   not  necessary  to   find  BB*.     We  get   the  required  results  by 
finding,    instead,    angles   a  and  c,    subtracting  these   from   w   to  get  angle 
b,    then  dividing  b  by  TT  . 

To  find   angle  a  we  will  make  use  of  the  fact  that  triangles  AND  and 

CBD  are  spherical  right   triangles.      Furthermore,    side   ¥   of  AND   is  — 9, 

Thus    if  we  can  find  side   A    then  Napier's  Rules   for  spherical   right  tri- 
angles will  give  us  angle  a.      A  similar  approach  is  used  to   find  angle  c. 
We  first   review  Napier's  Rules.      Given  a  spherical   right  triangle  as 
shown  in  the  sketch  below,   where  01 ,  rt  ,   and  tf    represent  the  central 
angles   subtended  by  the  respective  sides  of  the  triangle,   then  the  rules 
can  be  summarized  as   follows: 

Take  the  five  parts,   excluding  the  right  angle,    and  con- 
sider them  to  be  in  the  circular  arrangement   shown.      If 
"co-"  means    "complement  of"  then: 

(1)  the   sine  of   the  middle  part   equals   the  product  of 
the  tangents  of  the  adjacent   parts; 

(2)  the   sine  of  the  middle  part  equals  the  product  of 
the  cosine  of  the  opposite  parts. 
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To   find  angle  a  we  first   note  that  g   is   a   right   angle  since  A   is 
the  northernmost   point  of  the  orbit,    and  h   is   a   right   angle  by  con- 
struction.     The  sweep  width   is   2r   and  the  central   angle  subtended   by  r 
is  than  - 

From  Napier's   Rules, 

sin  £  -    cos(co-S)  cos(co-e) 


and  therefore 


Similarly, 


•     c  sin£ 

S'"S    B    IT^T    ■  (II) 


sin?  . 

SfaA    s     -7TZ7    '  Vz) 


sm  e. 

But 

X  +  S  »'£-*  ,  <i3> 

From(U),  (12),   and  (13), 

sin  V  sin S             sinV  cos  (*  +  \) 
Sin  A   -    r— ■ 7 -  <y-      r •  ( *+) 

After  some  manipulation  of  (14),  the  expression  for  X  is 

A  -  C°f    L        ^HpT)  sin? J   '     (  S) 

Again  using  Napier's  Rule, 

cos  a.  -  +**  ?  co+  \    .  (/*) 

From  ((5)  and  (16), 

-/f;    _-  r  *!*  z  t  gg  ^*  A) 5/"  >°7  7  (n) 
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To  find  angle  c  consider  the  following  sketch, 


Due  to  symmetry  between  northern  and  southern  hemispheres,  B*G  equals 
B* 'G' ' ;  therefore  angle  c  equals  angle  k.   But  angle  k  obviously  equals 
angle  f;  therefore  to  find  c  we  need  only  find  f .   Now  f  can  be  found  in 
the  same  way  we  found  angle  a  by  merely  using  spherical  right  triangles 
AND'1  and  C^'B^D' '  instead  of  triangles  AND  and  CBD.  Therefore  we  can 
use  equation  (17)  directly  by  replacing  (  A  +  §  )  with  (  A'  +  §').   But 
(  A'  +  £')  must  equal  Tf  -(  A  +  %    )  because  (  §  *  +  A  ')  equals  (  7T  -  •)   ) 
and  "i)  equals  (  X  +  §  )  by  symmetry.   Therefore,  after  the  appropriate 
changes  in  (17) ,  we  get 

sin  £  +  cos  (rr-  §  -A)  sin  fl  7      (is) 


C  =    cos 


-'[sec  r[- 


sin  (tr-S-X) 
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We  are  now  ready  to  obtain  an  expression  for  P(-d-  ,  <*  ,  r).   We 

said  P  would  be  equal  to  ,  where  b  is  obtained  from 

77 


b  -    7T-  (a  +  c)  .  (/<?; 


Using  the  relationships 

an  J    P(*>*>')=   £ 


in  equations    (17),    (18),    and   (19),   we  obtain      / 


/>(*,  «,r)=    J-  £fe**lj=x(*i'  *+  si»«  coi*>  + 


Cos 


-,rcsc-e;/s;n£_  s;nt<  cos  &)ji  (zo) 


On  a  reexamination  of  the  derivation  we  realize  that  (20)  is  valid  only 
for   <>(.£  &  -  %    .   Because  angle  a  is  zero  for  of.  >  &    -  £    ,    (20)  must 
be  modified  by  dropping  the  first  term  in  £*J  for  •£-£-<<£  ^-^S,  + 
Finally  P  must  be  zero  for  o( >  &  +  S    •   Therefore  the  expression  for 
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P  over  the  range     Q  £  u  £  ^£        is 

COS   [cos.*( 


P(W)=< 


IT  I  LCO&*\  /JJ     #-£*«<&+$ 


TT 


L 


O     ;        ++£*«*-£    . 


TT 


The  behavior  of  P  over  the  range  Oi«<-T-        is    illustrated   in 
figure  i. 
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APPENDIX   II 
Derivation  of  Estimate   3  for   Maan  Time  to   Next   Sighting 
Consider  the   following  sketch: 

N 


Wa  will  consider  our  line  of  sight  to  be  always  across  the  plane 
of  satellite  rotation,  and  let  the  earth  turn  relative  to  that  plane. 
I^t  D  =  time  for  single  revolution  of  earth  in  hours 

D  =  time  in  hours  for  earth  to  turn  far  enough  so  that  A  (fixed 
on  earth)  moves  twice  the  distance  QR.   (Thus  D'  is  the 
time  between  a  possible  sighting  on  an  upward  pass,  and  a 
possible  sighting  on  a  downward  pass.) 
TQ  =  Orbital  period  in  hours 
From  the  sketch, 

D'  =  arc  QR  D  =   angle  g  D 
arc  PR  tt 

But,  by  Napier's  Rule, 
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cos  g  =  cot  ■$■    tan  et    . 
Therefore  D'=     cos"1    (cot  »    tan  ^   )        ^ 

tr 

We  want  to  find  an  estimate  for  the  mean  time  to  next  sighting  of 
location  A,  given  A  has  suddenly  become  of  interest.   We  are  unable,  with- 
out great  difficulty,  to  describe  the  sequence  of  actual  sighting  times. 
Therefore  we  would  like  to  define  a  sequence  of  times  which  we  can  use  to 
approximate  the  sequence  of  actual  sightings,  and  upon  which  we  can  base  a 
reasonable  estimate  for  M(T^) . 

We  will  assume  that  the  sequence  selected  should  satisfy  two  condi- 
tions. 

Condition  1.   The  sequence  must  satisfy  the  condition  that  the 

mean  number  of  orbits  between  positive  orbits  in  n,  where  n  =  . 

P 

Condition  2.  The  sequence  should  be  proportional  to  the  sequence 
of  possible  sighting  times. 

It  is  obvious  why  condition  1  is  desired.  Condition  2  is  specified  be- 
cause it  will  require  the  sequence  of  assumed  actual  sightings  to  be 
spaced  in  a  way  which  suits  our  intuition.   Thus,  if  the  possible  sight- 
ings occur  in  a  sequence  such  that  the  time  between  possible  sightings 
has  the  sequence 

1  hr,  23  hr,  1  hr,  23  hr, , 

then  we  would  intuitively  expect  actual  sightings  to  be  spaced  such  that 
usually  t*o  sightings  would  occur  relatively  close  together,  followed  by 
a  long  period  with  no  sighting,  then  two  more  sightings  relatively  close 
together,  etc. 

Therefore  if  the  times  between  possible  sightings  has  the  sequence 
Z  =  D'  ,  D-D',  D'  ,  D-D'  , 
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then  we  want  a  constant  K  such  that  the  sequence  KZ  satisfies  our  two  con- 
ditions.  There  is  obviously  only  one  sequence  which  satisfies  these  two 
conditions  simultaneously.    Therefore  our  next  task  is  to  find  that  se- 
quence. 

To  find  the  sequence  KZ  we  need  only  find  a  K  which,  when  multiply- 
ing Z,  satisfies  condition  1,  because  any  sequence  KZ  where  K  is  constant 
will  automatically  satisfy  condition  2.   Now  the  sequence  Z  "cycles"  on 
every  second  sighting  possibility  with  a  cycle  time  of  D  and  therefore 
there  are  two  possible  sightings  in  a  cycle  of  Z.   We  would  like  KZ  to 
"cycle"  on  every  second  sighting  also,  but  the  cycle  time  must  be  equal  to 

2nT  ,  since  there  are  to  be  2n  orbits  per  every  two  sightings.   Therefore 

2nTn 
to  get  KZ  we  must  multiply  Z  by  —  .   Thus: 


fc  =.  £  n  %" 


We  then  have 


KZ=  {%£*',  if-^-n^:^-^.]. 


then, 


To  show  that  this  sequence  satisfies  condition  1  we  consider  a  single 

a* 
o 


cycle         (_     *    °  £>'    /?/7 /£ //? -D') )  The  mean  time  between   sightings    is 


//2»£j>'+  i?2L  (p-poj  =  nZ, 


2Ld 

and  therefore  the  mean  number  of  orbits  between  positive  orbits  is  n. 


We  have  fixed  the  mean  time  between  points  in  the  sequence  and 


the  "cycle"  of  the  sequence, 


36 


Having  defined  KZ,   the  sequence  of  times  between  assumed  actual 

it 

sightings,  we  can  say  that  assumed  positive  orbits  occur  at  times  T  , 

*   * 
T2,  T3,.... 


where 


r    -r* 


^*  = 


71*  *  **(f)         ;    i  odd 


ie.   [rj/J^T*,..-]  = 


{Z»T.@),  Z»T0,  **M?0  **,  -J 


To  estimate  M(T_  )  on  this  basis  consider  the  following  diagram: 
■ ZnX H 


~5- 


- — S. — - 


^-^-^T—Tsz  ~ 


T*  7T* 


where  T-  indicates  the  time  of  the  jth  assumed  sighting  and  the  time  be- 

"k  -k 

tween  T^  and  T-+2  is  defined  to  be  one  cycle.   Then  assuming  the  time  of 
occurrence  of  an  event  of  interest  at  A  is  uniform  over  the  satellite's 
life,  we  can  say  that  the  event  of  interest  is  equally  likely  to  occur 
at  any  time  T  contained  in  (T^,  T-+2)  • 

We  are  trying  to  estimate  the  mean  time  to  next  sighting  given  an 
event  of  interest  has  occurred  at  A.   However,  since  each  cycle  is  iden- 
tical, we  can  equivalently  find  the  mean  time  to  next  sighting  given 
T  €  (T*,  T. +2) •  The  density  of  T,  where  T  is  the  random  variable  denoting 
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.*      _* 


the  time  of  occurrence  of  the  event  of    interest   given  T  €   (T. ,   T.+~), 
is   then 


T  2  nT0 


Therefore 


*  /_.  \  Tq.     [mean  hme  +rom   T  +o   T^j  l 


7gz      Fmeun  Hnoe  from   T~  Jo  ^>^ 

ft,  [r*     ij*  +  n*l  ,    7e*  Ft*  _    £■*  +  %*' 7 

•k 

Substitution  of  the  appropriate  expressions   for  the  T-'s    in   (22)    results 
in 


and,    upon  substitution  for  D'    from  (21), 


[v  -  <»*-//£*+*  +*n  «)J2J .        (23) 


Thus,   we  obtain  a  general  estimate  for  M(T-r)    as   a  function  of   •&  ,  u    ,   and 
r. 
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It  is  of  interest  to  note  that  when  o<  =  0,  (23)  reduces  to 

M(TT)  = 


1/    2P   ' 
and  when  o(   =^,  (23)  reduces  to 

M(TT  )  =  — £-  . 
L     P 

These  limiting  forms  are  identical  to  estimates  1  and  2  respectively, 

which  is  a  desirable  result,  as  estimates  1  and  2  seem  appealing  at  these 

limits. 
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APPENDIX    III 
Comparison  of  Theoretical    and  Numerical    Estimates   for  MCT^) 
Our   analytical   results   are  compared  with   simulated   results    for  the 
following  data:      9   =   50°;    r   =   150n.nu;    TQ  =    1.6  hrs;    105  orbits. 

Q<  Analytical    Results  Simulated   Results 

0  32   hrs  35   hrs 

20  22.7   hrs  2  2.5   hrs 

40  16.5   hrs  16.7   hrs 

50  12.5   hrs  10.5   hrs 

The  simulated   results  were  obtained   from  a   Stanford  Research   Insti- 
tute study  which  computed  the   times   of   intersection   for   each  of   some   1800 
fixed   locations.    [2]      To  obtain  M(T,)    from  these  times   of   intersection 
for  each  location  the  occurrence  of   an  event  was   assumed  to   be  uniformly 
distributed  over  time  during  the   105  orbits   considered.      M(T-)    for   loca- 
tion  j  was   then  found  by  the  formula 

where       of  •  '}  L-f,2>''>  fir      are  the  times  of   intersection,    Si     is   the  total 
time  considered   for   location   J(  /~j   =  time  between  the   first   and   last    in- 
tersections  for    j)  .      This   formula   is   easily  verified   as   a   formula   for  the 
time  to    sighting  of   an  event   which  occurs   uniformly   in   (0,    /^j)   when  times 
of   sighting  occur  at   times    oft ,  o(2  ,    •••  0(n,  .     Then  the  overall  value  of 

M(T,  )   was   found  by  averaging  the  M(T-).      Since   locations  were  scattered 
h  3 

throughout  various    latitudes   only  those  within  a  small   band  were  considered 
to   find  M(Tt)    for  that  band.      For  example,    all   locations   between   39°   and 
41°  were  used  to   find  M(T  )    for  40°    latitude.      In  the  case  where  two   cities 

Li 
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were  so  close  together  that  considering  both  would  create  bias,    (Los 
Angeles   and  Long  Beach)    only  one  was   used. 

The  results  are  extremely  close  with  the  notable  exception  of  the 
range  where  we  would  expect  our  estimate  to  yield  the  best  results,    i.e. 
at     of  =  0  =  50°.      The  difference  here   is  most   likely  due  to  the  fact 
that   the  model    in  this  paper  assumes  the  satellite  sensor   "looks  straight 
down"  whereas  the  computer  model  has  the  sensor's  field  of  view  skewed 
toward  the  north  side  of  the  path.      This  difference  between  our  model 
and  the  computer  model   can  be  seen  from  the  following  sketch. 


symmetric    sensor 


/ 


s/ceweJ  sensor 


The  satellite  in  the  latter  case,  when  it  is  over  latitude  50°, 
could  have  its  sensor  centered  on  a  latitude  as  high  as  54°.  While  a 
four  degree  error  is  negligible  at  moderate  latitudes,  as  o<    approaches 
©  the  difference  becomes  quite  large.   If  the  fraction  of  positive  orbits 
from  the  computer  simulation  is  computed  and  used  in  our  formula  instead 
of  P(0,  o<  ,  r),  then  M(TL)  for  <*  =  50°  is  found  to  be  10.7  hrs.   There- 
fore the  variation  between  our  results  and  those  of  the  computer  simulation 
appears  to  be  caused  by  the  difference  in  fraction  of  positive  orbits  due 
to  slightly  different  physical  models. 
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APPENDIX  IV 
Derivation  of  Mean  Distance  to  Cross  Transmission  Area 

We  want  to  find  the  mean  distance  X  that  a  satellite  travels  when 
it  transits  a  transmission  area  of  a  ground  readout  station.  We  can 
estimate  this  distance  by  assuming  the  transmission  area  is  a  plane  cir- 
cle, and  that  the  track  of  the  satellite  is  a  chord  through  the  circle 
.with  the  distance  from  the  chord  to  the  center  being  uniformly  distribu- 
ted on  (0,r).   To  facilitate  the  discussion,  consider  the  following 
sketch. 


In  terms  of  the  sketch  we  want  to  find  E(X),  the  mean  of  X,  when 
the  random  variable  Y  is  uniform  on  AC,  that  is, 


From  the  sketch  we  get 


fc 


X'=    (r'-VT, 
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and,  because  Y  is  uniformly  distributed,  the  random  variable  X1  has 
the  distribution  function: 


.,.(£=£*)*■    K«r. 


(Zf) 


But  X  =  2x' ;  therefore  the  distribution  function  of  the  random  variable 


It  is  then  easily  shown  that 


TTr  (&) 


Equation  (26)  was  obtained  assuming  a  plane  circle,  whereas  the  actual 
transmission  area  is  a  segment  of  a  sphere.  The  result  is  a  reasonable 
approximation,  however,  because  the  transmission  areas  are  small  enough 
that  the  curvature  of  the  earth  is  negligible. 
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